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Introduction: Details of the differentiation of a 
global-scale lunar magma ocean (LMO) remain enig- 
matic, as the Moon is not simply composed of highlands 
anorthosite and a suite of mare basalts as inferred from 
early studies. Results from recent orbital missions, and 
the increasingly detailed study of lunar samples, have 
revealed a much larger range of lithologies, from rela- 
tively MgO-rich and “purest anorthosite” discovered on 
the lunar far side by the M? instrument on Chandraayan- 
1 [1] to more exotic lithologies such as Si-rich domes 
[2] and spinel-rich clasts [3,4] distributed globally. To 
understand this increasingly complex geology, we must 
understand the initial formation and evolution of the 
LMO, and the composition of the cumulates this differ- 
entiation could have produced. Several attempts at mod- 
elling such a crystallization sequence have been made 
[5-7], and have raised as many questions as they have 
answered. 

We present results from our ongoing experimental 
simulations of magma ocean crystallization, investigat- 
ing two end-member bulk compositions (TWM [8] and 
LPUM [9]) under fully fractional crystallization condi- 
tions. These simulations represent melting of the entire 
silicate portion of the Moon, as an end-member starting 
point from which to begin assessing the evolution of the 
lunar interior and formation of the lunar crust. 

Experimental: Starting compositions were synthe- 
sized from mechanical mixtures of anhydrous oxides. 
Experiments were performed in Depths of the Earth 
QuickPress piston cylinders at NASA Johnson Space 
Center, under nominally anhydrous conditions using 
graphite capsules and BaCO3 as a pressure medium. As 
this study aims to simulate crystallization from a magma 
ocean by cooling, experimental charges were heated to 
a superliquidus temperature before dropping to the tem- 
perature of interest. Experiments were held at tempera- 
ture for up to 72, and quenched by cutting the power to 
the assembly. Run products were analyzed by electron 
microprobe at 15kV and 20nA using a 10 um defocused 
beam on the glass phase. 

In order to simulate fractional crystallization, exper- 
iments with approximately 10% crystallization were tar- 
geted. The composition of the glass from each step was 
used as the starting composition for the next set of ex- 
periments, simulating 10% fractionation of crystals 
from melt. At each stage the total volume crystallized 
was calculated, and this value was used to calculate the 
pressure appropriate for the next set of experiments. 
Pressures below 0.5 GPa cannot be achieved in the JSC 


piston-cylinders, and so low pressure experiments were 
run in Deltech CO-CO2 gas mixing furnace buffered at 
IW. 

Results: The full fractional crystallization sequence 
is complete in the LPUM composition, and 74% com- 
plete in the TWM composition. A weighted average of 
the cumulate pile composition based on modal abun- 
dance of each phase, taking into account the changing 
composition of a given mineral as fractionation pro- 
gresses, has been calculated, and is within 5% of the 
LPUM starting composition, indicating mass balance 
throughout the experimental sequence. The experimen- 
tally derived LMO (Fig. 1) is characterized by the pre- 
domination of olivine early in the crystallization se- 
quence. Olivine is the only crystallizing phase for the 
first 50% of LMO crystallization in both bulk composi- 
tions. 
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Figure 1: Experimental results of fractional crystallization 
experiments on two LMO bulk compositions. 

After 50% crystallization, the liquidus phase 
switches to low-Ca pyroxene in the LPUM composition, 
whereas both low-Ca pyroxene and Cr-rich spinel co- 
crystallize with olivine in the TWM composition. At 
75% crystallization plagioclase becomes a liquidus 
phase in the LPUM composition, but at 74% crystalli- 
zation in TWM (the highest degree of crystallinity yet 
explored), plagioclase is not present, despite having 
more Ca and Al in the bulk composition than does 
LPUM. The LPUM cumulate pile consists of approxi- 
mately 10.5% modal plagioclase, which corresponds to 
a crustal thickness of approximately 60km. 

Discussion: The experimental results indicate that 
neither end-member bulk composition will produce a 


cumulate pile with plagioclase abundance that is propor- 
tional to a crustal thickness of 45 km or less. LPUM pro- 
duces plagioclase abundance consistent with a crust of 
~ 60km, and although experiments on TWM are not 
complete, we expect it to produce a thicker layer of an- 
orthosite due to the enrichment of refractory elements in 
this composition — the higher abundance of Ca and Al 
were initially proposed to account for a crust of up to 
100 km depth. However, our model does not take into 
account trapped interstitial liquid. [6] showed that in- 
creasing the amount of trapped interstitial liquid in their 
model decreases the amount of anorthosite in the final 
cumulate pile as Al-rich liquid is sequestered in the 
mantle. Furthermore, the depth of the magma ocean and 
the mode of crystallization are important factors in the 
amount of plagioclase produced. A shallower magma 
ocean would likely stabilize plagioclase earlier in the se- 
quence, as it the major aluminous phase stable at low 
pressure. At higher pressures plagioclase crystallization 
is suppressed, and Al-rich pyroxenes, spinel and even- 
tually garnet are stabilized. There is evidence of this in 
the experimentally produced TWM cumulates, which 
contain a spinel-structure phase which is a solid solution 
between chromite (FeCr2O4) and spinel (MgAlI2Ox.) 
from 55 — 70% crystallization (~1.5 — 1 GPa). 

Alternative LMO crystallization models. As alluded 
to above, several groups advocate for a shallower 
magma ocean rather than full mantle melting [e.g. 10], 
despite evidence suggesting the heat flow involved in 
accretion would be more than enough to account for 
whole Moon melting (9,11]. Given the stability of pla- 
gioclase over other aluminous phases at lower pressure 
however, a shallower magma ocean would likely in- 
crease the thickness of the plagioclase crust produced. 

Another alternative model for LMO crystallization 
is that of a period of equilibrium crystallization, in 
which crystals remain entrained in the magma, followed 
by fractional crystallization after a critical melt viscos- 
ity is reached [12]. This two-step approach will produce 
a less extreme fractionation of the magma, and poten- 
tially somewhat less plagioclase, and therefore could be 
a preferred mechanism for LMO crystallization. Com- 
putational models by [5] have been the framework for 
LMO crystallization in the lunar community, but exper- 
imental work by [13] called some of the conclusions 
into question. 

Although experimental work on the fractional part 
of LMO crystallization following initial equilibrium 
crystallization by [13] are ongoing, the MELTS algo- 
rithms [14,15] allow us to investigate the likely crystal- 
lization sequence for a residual liquid after 50% equilib- 
rium crystallization. There are limitations to using 
MELTS, and there is a lack of experimental data to cal- 
ibrate the model for calculations at the P,T,x appropriate 


to LMO crystallization. Nevertheless it provides a use- 
ful tool to investigate crystallization sequences before 
embarking on a full experimental program. We per- 
formed fractional crystallization calculations on the 
LPUM 50% residual liquid from [13] at 2, 1, 0.5 GPa 
and 1-bar to simulate different depths of magma ocean. 
Plagioclase is present only in the calculations at < 1 
GPa, and forms approximately 10% of the cumulate pile 
— similar to our experimental outcomes. The next step is 
to carry out a polybaric fractionation sequence, as at 
higher pressures garnet is stabilized. The geochemistry 
of some lunar basalts suggests that garnet is residual in 
their source region [16]. Experimental work to date has 
failed to produce significant garnet in cumulates at an 
appropriate depth for basalt source regions, but crystal- 
lization under hybrid equilibrium/fractional conditions 
at higher pressure could allow for garnet crystallization, 
which would lead in turn to less plagioclase crystalliza- 
tion, consistent with a thinner anorthosite crust. Alt- 
hough further experimental and sample work is neces- 
sary, it is becoming increasingly clear that a simple 
model of LMO crystallization is no longer adequate to 
explain lunar geology. 
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Figure 2: MELTS model outcomes of fractional crystalli- 
zation of the residual liquid after 50% equilibrium crys- 
tallization of LPUM [13]. 
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